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Abstract 

With a consistent definition of transverse-momentum-dependent(TMD) light-cone wave func- 
tion of i3-meson, we show that the amplitude of the radiative leptonic B-decay can be factorized 
at one-loop level as a convolution with the wave function and a perturbative coefficient function, 
combined with a soft factor. In this TMD factorization, the transverse momenta of partons in 
the i?-meson are taken into account and all soft divergences are contained in the wave func- 
tion and the soft factor. The coefficient function is infrared-safe. The factorization works on a 
diagram-by-diagram basis and is possible to extend beyond one loop. With the factorization the 
large logarithms in the perturbative function can be simply resummed. Our work shows that 
the result of collinear factorization for the decay can be derived from that of TMD factorization. 
Therefore, the two factorizations for the case here are simply related to each other. 
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1. Introduction 



Exclusive B-decays play an important role for testing the standard model and seeking for new 
physics. Experimentally they are studied intensively. Theoretically, there are two approaches 
of QCD factorization for studying these decays. One is based on the collinear factorization pp, 
in which the transverse momenta of partons in a B-meson are integrated out and their effect 
at leading twist is neglected. The collinear factorization has been proposed for other exclusive 
processes for long time[2]. Another one is based on £;T-factorization[31 or pQCD approach, where 
one takes the transverse momenta of partons into account at leading twist by means of kx- 
dependent light-cone wave function. We will call such a factorization as transverse momentum 
dependent(TMD) factorization. The advantage of the TMD factorization is that it may eliminate 
end-point singularities in collinear factorization [1] and some higher- twist effects are included. The 
knowledge of the transverse momentum dependent (TMD) light-cone wave function will provide a 
3-dimensional picture of a .B-meson bound state. However, it was not clear how to define the TMD 
light-cone wave function in a consistent way to perform a TMD factorization because of light-cone 
singularities (Sj- 

Similar problems also appear in defining TMD parton distributions and fragmentation functions 
if one tries to do TMD factorization for inclusive processes. In general the light-cone singularities 
appear if a parton emits gluons carrying momenta which are vanishingly small in the +- direction 
but large in other directions in a light-cone coordinate system. In a collinear factorization for 
an exclusive or inclusive process, these singularities are cancelled between different contributions 
if the transverse momentum of the parton is integrated out. If the transverse momentum is not 
integrated, the singularities are not cancelled. 

For inclusive processes like Drell-Yan, semi-inclusive DIS etc., it has been shown that one can 
consistently define TMD parton distributions by using gauge links in the direction off the light- 
cone direction and the TMD factorization of inclusive processes can be done without light-cone 
singularities El El- The TMD parton distributions defined with these gauge links will depend 
on the deviation of the direction from the light-cone direction. The evolution with this depen- 
dence is controlled by the Collins-Soper equation 6 j which leads to the so-called CSS resummation 
formalism [7JIE1IH1- This formalism is for resummation of large logarithms appearing in the collinear 
factorization. In that sense the TMD- and collinear factorization are related to each other. But 
the similar relation in exclusive B-decays has not been studied. 

We have proposed in [10 a to consistently define the TMD light-cone wave function of B-meson 
by using gauge links off the light-cone direction and studied its relation to the usual light-cone 
wave function in the collinear factorization. With the consistent definition it is important to show 
that the TMD factorization can be consistently performed. The relation between two factorization 
approaches may be then established. As a first step towards these goals we study in this paper TMD 
factorization for the radiative leptonic decay of B-meson. We also plan to study TMD factorization 
for B — > 7r form factor and other decay processes. It should be noted that the definition of the 
TMD light-cone wave function is not unique, different definitions are possible. A different definition 
can be found in With different definitions the most important thing is to show that one can 
perform TMD factorization consistently with one of these definitions, at least at one-loop level. To 
our knowledge, there is so far no such a study beyond tree-level for exclusive B-decays. We will 
show that with the definition given in ^01 the factorization can be done at one-loop level for the 
process studied in this paper. 

The radiative leptonic .B-decay has been studied extensively |12l IT51 ITU ITdI [TBI ITS] . The effect 
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of strong interaction in the decay is parameterized with form factors. These form factors have been 
studied in ^2l El El El El El with QCD factorization. It has been shown that the form factors 
can be factorized as a convolution with a perturbative coefficient function and the light-cone wave 
function of the Z?-meson in the collinear factorization|14 | I15 | I16[ IT8] . In these works the transverse 
momentum of partons is integrated out. It results in the convolution only with the +- component 
of the parton momentum. In ^3] the transverse momentum of the parton is not integrated and 
is explicitly taken into account, but the consistency of the definition of the TMD light-cone wave 
function is not addressed and the problem of the gauge invariance of the definition is ignored. In 
|17j the decay is studied with /c^-factorization or TMD factorization, but the TMD light-cone wave 
function employed there has the light-cone singularity. 

With our gauge-invariant definition we can show with TMD factorization that the form factors 
take the factorized form: 

(/> + ®S®H. (1) 

In the above (/>+ is the TMD light-cone wave function, S is a soft factor, H is a coefficient function 
which can be calculated with perturbative QCD and is free from soft divergences. and S are well- 
defined matrix elements of QCD operators. The convolution here is not only with +- components 
but also transverse components of parton momenta. In this paper we prove the factorization at 
one-loop level. We show that the cancellation of all soft divergences is on a diagram-by-diagram 
basis. This is important for extending our factorization beyond one-loop level. In the case studied 
in the paper, the TMD factorization is similar to the collinear factorization because there is no 
parton or hadron in the final state. But it is important to show first that the TMD factorization 
works for this simple case and then extend the TMD factorization to other cases. An interesting 
fact with the TMD factorization is that it provides a simple way to resum large logarithms in the 
perturbative function H, as we will show in the paper. 

As mentioned before, TMD factorization for an inclusive process can be related to the cor- 
responding collinear factorization. One can expect that such a relation also exists for exclusive 
-B-decays. Indeed, in the case studied here, such a relation exists and it is simple: Both factor- 
izations are equivalent, i.e., one can derive the result of the collinear factorization from our TMD 
factorization. We will show this in this work. One reason for this simple relation is that there is 
no hadron, hence any parton in the final state. 

Our paper is organized as the following: In the next section we define the TMD light-cone wave 
function and give its one-loop result in detail with a general partonic state, which will be used to 
perform TMD factorization. A detailed discussion about the TMD light-cone wave function and 
its relation to the usual light-cone wave function in collinear factorization can be found in ^U]. In 
Sect. 3 we introduce our notation for the decay and the result of the factorization at tree-level. 
In Sect. 4 we will complete the factorization at one- loop level and determine the soft factor. In 
Sect. 5 we show that the result of the collinear factorization can be derived from that of the TMD 
factorization and establish the relation between the two factorizations for the decay. In Sect. 6 we 
will make an attempt to re-sum large logarithms in TMD factorization. Sect. 7 is our conclusion 
and outlook. 

2. A Consistent Definition of the TMD Light-Cone Wave-Function 

In this section we give our definition of the TMD light-cone wave function and its one-loop 
result in detail with a general partonic state. A brief report of the result and the study of the 
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relation to the light-cone wave function in the collinear factorization can be found in jlUj . 

We will use the light-cone coordinate system, in which a vector is expressed as a M = 
(a + ,a~,a ± ) = ((a + a 3 )/ v / 2, (a - a 3 )/\/2, a 1 , a 2 ) and a\ = (a 1 ) 2 + (a 2 ) 2 . For b-quark we will use 
the heavy quark effective theory (HQET). To define the TMD light-cone wave function we introduce 
a vector u M = (u + ,u~,0, 0) and the definition is given in the limit u + « u~ jl()j: 

Mk + ,kx,C,^)= J ^^e lfc+ ^-^^M0|^)Lt(oo,z) 7 + 75L u (oo,0)/ l (0)|S(t;))U +=0 , (2) 

where h(x) is the 6-quark field in HQET. and L u is the gauge link in the direction u: 

L u (oo, z) = Pexp ^— ig s J dXu ■ G(Xu + z)^j . (3) 

In the above, the 5-meson moves with the velocity = (v + , v~ , 0, 0), i.e., in the z-direction. The 
limit should be understood that we do not take the contributions proportional to any positive power 
of u + /u~ into account. This definition is gauge invariant in any non-singular gauge in which the 
gauge field is zero at infinite space-time. It has not the mentioned light-cone singularity as we will 
show through our one-loop result, but it has an extra dependence on the momentum k + through 
the variable C 2 = 4(u • k) 2 /u 2 = C{k + ) 2 , or an extra dependence on £ 2 . The evolution with the 
renormalization scale /i is simple: 

V q- = (7g+70)^+(* »*±»C»M)» ( 4 ) 

where 7 g and jq is the anomalous dimension of the light quark field q and the heavy quark field 
h in the axial gauge u ■ G = 0, respectively. In the remainder of the paper we will not indicate 
the //-dependence explicitly if it does not cause any confusion. It should be noted that one can 
not simply relate 0+(A; + , k±, £) by integrating k± to the light-cone wave function in the collinear 
factorization, whose definition can be found in [TJ||. The reason for this has been discussed in detail 
in P3J. 

To perform TMD factorization one needs to calculate the wave function with perturbative QCD, 
in which the i?-meson is replaced by a partonic state. We take the partonic state \b(mbv+kb),q(k q )) 
to replace the S-meson in the definition, the momenta are given as 

kq = (k q , k q , k q ±), kfr = (k b , k b , — k q ±). (5) 

These partons are on-shell, i.e., k 2 = m 2 and v ■ k b = in HQET. It should be noted that we take 
a finite k q ± without loosing generality. The quark mass m q will regularize collinear singularities. 
We also introduce a gluon mass A to regularize infrared singularities. The variable k + of the wave 
function is from to oo in the heavy quark limit. Actually, from the momentum conservation, it is 
from to P + = mi ) v + + k b + Under the limit mj, — > oo we have P + — > oo. As discussed in jlUj . 
if we set P + to be oo at the beginning, it may result in some ill-defined distributions. Therefore 
we should take a finite P + in the calculation and take the limit P + — > oo in the final result. For 
results obtained in this paper we will take the limit where it does not introduce any problem. 
At tree-level, the wave function reads: 

</>^(k + ,k ± ,0 = v(k q ) 7 + l5 u(k b )5(k + - k+)5 2 {k L - V). (6) 
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We will always write a quantity A as A = +A^ H , where A^ and A^ stand for tree-level- 

and one-loop contribution respectively. At one-loop one can divide the corrections into a real part 
and a virtual part. The real part comes from contributions of Feynman diagrams given in Fig.l. 
The virtual part comes from contributions of Feynman diagrams given in Fig. 2., these contributions 
are proportional to the tree-level result. 




(a) (b) (c) (d) 



Figure 1: Diagrams of one-loop contributions. Thick lines stand for 6-quark, double lines represent 
gauge links. 

To illustrate how to calculate these contributions and how the limit u + « u~ is taken, let 
us consider the contribution from Fig.lc. After integrating the — component of the momentum 
carried by the exchanged gluon the contribution reads: 

2a, 



i)+(k + ,k ± ,()\ic = -^v(k q )j + -f 5 u{k b )u-v6(k + - k+) 

2q+ 



v+{q\ + A 2 ) + 2v-(q+) 2 + iO u+{q\ + A 2 ) + 2u'{q+) 2 + iO' 



(7) 



with q + = fc+ — k + and q± = k q ± — k±. If we simply set u + = 0, the contribution is proportional 
to l/(k + — k^) and divergent at k + = k^ . This is the mentioned light-cone singularity. With the 
nonzero u + the contribution is finite for any k + . The contribution in the limit u + « u~ can be 
derived by taking the contribution as a distribution of k + and it reads: 



2a s _ 

4> + (k + ,k ± ,()\ic = -^v{k q )-i + ^u{k b ) 

6{k+ - k+) 



q+(ql + \2 + (2(q+)2) 



+ o(C 2 ), 



c 2 = 



2u 



u 



C 2 . 2r 

' (k+) 2 ' 



+ 



(8) 



In the above the limit P + —> oo is already taken. From the result we can see that the light-cone 
singularity is regularized with the finite but large £«• The other contributions of the real part are: 



(j) + {k + , fc_L,C)|la 
Fl 

4>+(k + , All, C)| 16 



2a 



-^v(k q )"f • v(7 • (q - kq) + mq)"f + <y 5 u(k b )F 1 
dq- 1 1 



1 



2tt (q - kg) 2 -m 2 + iO q 2 - A 2 + iO v-q + iO' 

k+ (e{k+-k+) 



^v(k q )~f + -f 5 u(k b ) 



kt-k+ 



+ 



1 



+ 



2{q\ + A 2 ) 



5(k~ 



c 2 



5 



A g = k+((q ± -xk q± ) 2 +x 2 m 2 q + (l-x)\ 2 ), x = 1 - — , 

2a 1 
<M& + ,fc-L,C)|id = -7^v(k q )j + j 5 u(k b )5(k + - £;+) ^-, + ^ 2 . (9) 

The integral i*i for the contribution from Fig. la can be done easily, but it results in a lengthy 
expression. We will show later that the contribution will not affect the perturbative coefficient 
function H. 











(a) 


(b) 


(c) 


m ^0 m 


"* — ^Cf 





(d) (e) (f) 



Figure 2: The virtual part of the one-loop correction. 



The virtual part of the one-loop correction is from the Feynman diagrams given in Fig. 2. 
Contributions from each diagrams are: 

,2 

Mk + ,k ± ,()\2a 



cf>+(k + , k±,Q\2b = Mk + , k ± , C)| 2e = <$\k + , k ± , C) — In ^, 



4>+(k + ,k ± ,()\2c 
Mk + ,k ± ,()\2f 
<P+(k + ,k ± ,C)\ 2d 



0l O) (fe + ,^,C)-^ln^ln 



^\k+,k ± ,0- 



a 



3vr A 2 ( 2 ' 

2 a2 



6tt 



TO 2 ^2 n„2 

2 In r- + 2 In - In 2 At - 2 In In C 

A 2 



c 2 



c 2 



rn 



^\k + ,k ± ,C)-^ 



The complete one-loop contribution tjyr} is the sum of contributions from the 10 Feynman diagrams 
in Fig.l. and Fig. 2. With these results one can derive the evolution of £. For this we transform 
the wave-function into the impact parameter 6-space: 



mi 



I1K, 



+ 4 



In At + 21n— % -4 



mi 



A 2 



(10) 



+ (fc + ,&,C,/z) = I d 2 k ± e ik ^ b Mk + ,k ± ,(,n), 



the evolution reads: 

d 



4a Sl (WH 2a Sl /x 2 e 

m In — 

3vr 4 3vr C 2 



(11) 



(12) 



The first factor is the famous factor K + GjHlE], ^ ne ^ as ^ factor comes because we used HQET for 
the heavy quark. 
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Before ending the section, we briefly discuss the heavy quark limit P + — * oo. For the usual 
light-cone wave function, this limit will result in that the wave function is not normalizable as found 
in |19l I2U| and it is shown through an explicit calculation with perturbative theory in ^0] . For the 
TMD light-cone wave function it is normalizable if the transverse momentum is not integrated out. 
When we transform the TMD light-cone wave function into 6-space, we should keep P + finite. 

3. Notations and Factorization at Tree Level 

We consider the radiative decay of the S-meson B which contains at least a 6-quark and a light 
anti-quark q: 

B^-i + t + U. (13) 

We take a frame in which B moves in the z-direction with the velocity = (v + ,v~,0, 0) and the 
photon with the momentum p^ = (0,p~, 0, 0). It is worth to mention here that this decay has not 
been observed so far. An upper bound for the branching ratio is given in |21| : 

Br(B -^j + £ + u) < 2.0 x 1(T 6 . (14) 

In the decay the effect of the strong interaction is controlled by a matrix element of the operator 
<77 M (1 — 75)6 with b(x) being the 6-quark field in the full QCD. Since we use HQET for the heavy 
6-quark, the matrix element can be matched to HQET: 

(l(e*,p)\q^(l ~ 7S)6|B(«)) = f(v)(l(e*,p)\q^(l - ^ 5 )h\B(v)}, (15) 

where f{n) is the matching coefficient. It is given by: 

/(M) = l + ^(31n^-2)+0( o J). (16) 

The HQET matrix element can be parameterized as 

T» = -^(7{e%p)\<tf(l-^)h\B[v)) 
V47ra 

= e^e* u v pPa F v (v ■ p) + % (v ■ pt? - v ■ eV) F A (v ■ p). (17) 

In the rest frame of B, v -p is the energy of the photon. The invariant v -p can be from to Mb/2. 
The photon is emitted by quarks inside the S-meson. If v ■ p is large, i.e., v ■ p 3> Aqcd those 
quarks will change their momenta significantly, i.e., the emission becomes a short-distance process. 
This leads to that those form factors, hence the matrix element can be studied with perturbative 
QCD, in which one can separate short-distance- and long distance effect by factorization. 

To show the factorization, one usually replaces hadronic states with reasonable parton states, 
then calculate processes which need to be factorized and nonperturbative objects like wave functions 
in our case to extract the perturbative coefficient functions. A factorization means at least that 
those coefficient functions do not contain any soft divergence. For our purpose we replace the B 
state \B) with the partonic state \qb). The momenta of the partons are the same as in Eq.(5). At 
tree- level the contribution to the matrix element is given by the two diagrams in Fig. 3. The second 
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Figure 3: Tree-level contribution to the matrix element. The thick line is for the 6-quark, the black 
dot denotes the insertion of the operator. 



diagram will not contribute in the heavy quark limit by noting the fact v ■ e* = for a real photon. 
The tree- level amplitude T M reads: 

?f = QMhh • ■ " 75)n(fc&) ' (18) 

where Q q is the charge fraction of q. In TMD factorization one will neglect the transverse mo- 
mentum of initial partons in nominators of propagators but keep it in the denominators. The 
case studied here is rather special because the denominator does not depend on the transverse 
momentum. With a little algebra one can show that 



iQ 



2v ■ p 



e^ vpa e w v p p a + i (v • pe* - v ■ e*p^ ■ -^v(k q )>y + -y 5 u(k b ) H , (19) 



F V = F A = J dk+d 2 k ± Mk + ,k ± ,C,^- (20) 



where • • • denotes the neglected /^-dependence from the quark propagator and the contribution 
from the partonic state which does not have the same quantum numbers as B does. 

With the tree-level result of the TMD light-cone wave function, we obtain the factorization for 
those form factors: 

2v ■ p . 

At the orders considered in the work, Fy is always the same as Fa- We will write our factorization 
formulas as: 

F V = F A = I dk + d 2 k ± dl + d 2 l ± Mk + , k ± ,()S(l + , l±, ( u )0(k + + l + )H{k + + l+, Q, (21) 

zv ' p j 

so that at the leading order of a s the perturbative coefficient function H and the soft factor in 
perturbation theory at tree-level reads: 

# (0) (^ + ,C«) = S^(k + ,k ± ,( u ) = 5(k + )5 2 (k ± ). (22) 

It is noted that at the leading order H does not depend on k± in the case studied here, while in the 
other cases like B — > ir transition it does. If one replaces the B-meson state with a partonic state 
of off-shell partons, one can have a which depends on £^|T7]- But the amplitude with the 
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state of off-shell partons is not gauge-invariant. In general it is not clear if the factorization with 
such a state can be made in a gauge invariant way. 

At tree-level one can not determine the form of the soft factor S, because it is designed to 
subtract infrared divergences at higher orders of a s . It should be a 5- function at tree level. At 
one-loop level with the partonic state, the factorization formula takes the form: 

T« ~ cpf ® 5< > ® # (1) + <^ 0) ® ® + $ ® §W ® (23) 

the soft factor should be chosen so that all soft divergences of Tfp are contained in the second- 
and third term and is free from any soft divergence. The soft factor should also be chosen so 
that the factorization can be extended beyond one-loop level. 

4. The Soft Factor and Factorization at One-Loop Level 

In this section we will perform TMD factorization at one-loop level and determine the operator 
form of the soft factor. The perturbative coefficient function will also be determined at one-loop 
level. Let us first consider the one-loop corrections to the amplitude T^. The corrections are from 
diagrams given in Fig. 4. 




Figure 4: One-loop contribution to the matrix element. The thick line is for the 6-quark, the black 
dot denotes the insertion of the operator. 



The contribution from Fig. 4a reads: 

T^Ua = iQ g g s C F j7m v y k in ■ v T — — tt~^' € 



(2vr) 4 v q ' ' {k q + If -m 2 + is' (p - k q - I) 2 -m 2 q +ie 



•7 M (1 - 75) -u(k b ) 

-V ■ I + IE I 2 - X 2 + IE 

.■in I dH -11 ^ . + Q 



' ' ' +0(E~ 2 ). (24) 



-2p ■ (k q + l)+i£ -v-l + ie I 2 - A 2 + ie 

It is easy to find that this contribution up to a power correction is exactly represented by the 
contribution from </>+|i a to the third term in Eq.(23). Therefore, this contribution and 0+|i a is 
irrelevant for the determination of and . Also, the contributions from Fig.4e and Fig.4f 
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to are reproduced by the contributions from 



respectively. The contributions from other diagrams to T^> are: 

2k Q -p~ 



1 2a and 4>+\2d in the third term in Eq.(23), 



■lb 



T (0) . ^£ 

" 3vr 



+ 21n- 



1 pL\ic 



Id 




-7T 



r(°) • 
" 3vr 



(25) 



Our T 1*1 4b agrees with that of (14] . but is not in agreement with that in ^3]. The other contributions 
are in agreement with In these contributions there are no infrared singularities. They have 
only a collinear singularity from Fig. 4b, represented by lnm,j. The relevant contributions to <fi^ (g) 
5(0) ® H (0) by using are: 



dk + d 2 k±— [4> + \ lb (k + , k±) + (p + \ 2e (k + , k±) + ^ + | 2 /(A; + , fcjj] /v(k q )-/ + j 5 u(k b ) 



37rA; 



2 + In 



In 2 



c 2 

+ 21n^ +ln 



A 2 



-7T 



d 2 £; ± 



1 



7T A 2 + k\ 



In 



A 2 + fc 2 



c 2 



Wr 



J dk + d 2 k L -^ [(f>+\ic{k + , fej.) + 0+| 2 6(fc + , fe±) + 0+|2c(A: + , fcj.)] M^h + 75? 







37T/Cq" 



ln >7 + ln >7 ln ^| 
i 



-7T 



In 2 



d 2 ki 



1 



d/c + (i 2 fc_|_^-0 + | 1( i(A; + , k±)/v(k q )~/ + j 5 u(kb) 



2a s 
3irka 



7T A 2 + £: 2 
c2 2 A; 1 



In 



A 2 + fc 2 

c 2 



7T A 2 + A: 2 



(26) 



Comparing the above two equations, we find that the collinear singularity from Fig. 4b. is re- 
produced by the contribution in W from Fig.2f. But, there are many infrared singularities in 
cffl <g> £K°) ® F(°) ~ W a + W& + W c + W rf + W e , where W a is the contribution from Fig.la and 
W e is the sum of contributions from Fig. 2a and Fig. 2d. There are even ultraviolet divergences. 
However these divergences are closely related to corresponding infrared singularities, as they stand. 
Once these infrared singularities are subtracted, one can expect that those ultraviolet divergences 
are subtracted too. As mentioned before, the contributions of W a and W e represent those to the 

from Fig. 4a, Fig.4e and Fig.4f. To complete the factorization one needs to find the soft factor 
so that all infrared singularities in W b + W c + and also the divergent integrals over k± are 
subtracted by the soft factor. 

Clearly all these infrared singularities are from the TMD wave functions, i.e., from contributions 
from Fig.l. and Fig. 2. By using the eikonal approximation one easily finds that these singularities 
can be reproduced by the expectation value of the product of gauge links: 



5 4 (z~,z_l) 
Lu(z, -oo) 



dz-d z z ± e iq z - iq± ' z± S 4 {z-,z ± ), 
^Tr(0|r[4(^,-oo)4(oo^)L u (oo,0)L„(0,-oo)] \Q)\ z+=0 , 
Pexp (^—ig s J d\u ■ G(Xu + z) \ , 
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L v (z, — oo) = Pexp^—ig s J dXv ■ G(Xv + z)^j , 



(27) 



where the direction in L„ is chosen as u + » u~ . The gauge link L~ just simulates an anti-quark 
q in the initial state and L v the 6-quark in the initial state. If one only takes this product of the 
gauge links into account, one can expect that the quantity 



S±(z~, b) 



(28) 



is free from infrared singularities. This is checked at one-loop level. However, because part of 
contributions from </>+\ which are from Fig. la, Fig. 2a and Fig. 2d, is already used up to subtract 

soft divergences in Tjp , as discussed before, one can not expect that the soft factor S can be formed 
with S4 only. We will turn to this point later and concentrate at moment on the perturbative results 
for S 4 . 

At tree-level, the result is just a (^-function: 



(29) 



At one-loop level, there are contributions from diagrams which have a one-to-one correspondence 
to those diagrams given in Fig.l. and Fig. 2., in which one only needs to replace the light-quark 
line with the double line of the gauge link L~. The corresponding contributions as a distribution 
of q + for the range — k + < q + < 00 under the limits u + — ► and u~ — ► are: 



£4 (g + , <?_]_) I la 

S4(q + ,q±)\ib 
£4 (<? + , <?_!_) |lc 
S4(q + ,q±)\id 



2a s 
3tt 2 

2a s 
3tt 2 



0(q 



q + {ql + \ 2 + Q{q + ) 2 ) 
1 



0(- 



ql + V + Cl(q + ) 2 
2a, / %+) 
3vr 2 [\q+(ql + \ 2 + C 2 (q + ) 2 ) 
2a s 6{q+) 



q + (qi + \ 2 + (l(q + ) 2 ) 



1 



+ " imaginary part" , 

c 2 



In 



+ 



+ 2<% H 



1 



'ql + X 2 ln C 2 



(?i + A 2 g 2 +A 2 

c 2 



3vr 2 g 2 +A 2 ' 

and the contributions from the diagrams corresponding to those in Fig.2 are: 



(30) 



S4(q + ,q±)\2a 

Si(q + ,q±)\2c = Sr(q + ,q±) 
S4(q + ,q±)\2f = S!?\q + ,q ± ) 



S4(q + ,q±)\ 2 b = S 4 (q + ,q ± )\ 2 d = S 4 (q + ,q±)\ 2 e = Sf ] (q + ,q±) 

2 a2 



3^' 



:ir: ~ A 2 C 2 ' 
2 rl 
In ^ In ■ , 



3tt A 2 



with £~ = 2u /u + . If we identify the soft factor S(z ,b) as S\ 



0^ (8) £W (8) F (0) can be grouped similarly as those to <g> <g> They are 



(i) 



(31) 



b), their contributions to 



U a 



J 



dl + d 2 l 1 



4 + /+ 



Sl\la(l + ,l±) 
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-In 



3irkq I 2 



dl + d 2 l± 



A 2 
1 



57T 2 
IT 



-In 



kt + V 



37^ 



1 9 

-A + - In 2 — j 



2"" A 2 

[S^llb + S^be + 5*4|2/] (/ + , ?1 
\ 2 ii 2 ii 2 t 2 

■ ^-ln^-ln^ 
A 2 m A 2in C 2 



+ A > + " imaginary part" , 



In 



d 2 k ± 



1 



vr k 2 , + A 2 



In 



A: 2 +A 2 



dl + d 2 l ± 



1 



37r/co 6 



• -' • - hi 2 



[Si\lc + S4I26 + £4 1 2c] (Z + , Zj_) 



x 2 n 2 1, 2 r 2 

_ ln _ _ ln _ lri _ + 



1 



vr /c 2 + A 2 



In 



fc 2 +A 2 ' 
C 2 



17« 



dl + d 2 l± 



1 



cZZ+d 2 Z 



1 



;^4|idG + ,/±) 



2q s 
37tA:o 



d 2 /fc 



1 



vr /c 2 +A 2 ' 

2a s /i 2 
~~37r~ A 2 ' 



where A is a divergent quantity: 

A = 



&+ /■ ci 2 fc ± 2 
k ^o Ln k+J vr C 2 (A;+) 2 + fe 2 ' 



(32) 



(33) 



which will be cancelled in U a + Comparing the sum U a + Ub + U c + Ud + U e with Wb + W c + W<2, 
we note first that the divergent integrals over k± in Wb, W c and Wd are completely subtracted by 
those in Ub, U c and Ud, respectively. Also the infrared singularities with In A in Wb, W c and Wd are 
completely subtracted by those in Ub, U c and Ud, respectively. The remaining infrared singularities 
are only from U a and U e . 






Figure 5: One- loop contribution to S2. The double lines represent the two gauge links. One is for 
L v , the other one is for L~. 



These remaining singularities can be reproduced by the product of the gauge links: 

S 2 = ^Tr(0|T [4(0, -oo)L v (0, -00)] |0). 



(34) 



At leading order = 1. At one-loop level, the contributions are from the diagrams given in 
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Fig.5. 



'2 1 5a 



^ln^ln s < 
3vr A 2 



c 2 

_ 

c 2 



+ " imaginary part" , 



S 



2\5b 



S 



2 1 5c 



37 ln A^- 



(35) 



Now we turn to the imaginary or absorptive part. In the amplitude it has an absorptive part 
from the box diagram Fig. 4a and its contribution is already contained in the contribution from the 
TMD wave function in Fig. la. The remaining parts can not have an absorptive part. Therefore, 
one should eliminate possible absorptive part in the soft factor. At one-loop level, the imaginary 
part from £4 is the same as that from Si. But this is from perturbative theory. To eliminate the 
absorptive part one can simply take the real parts of those products of gauge links. Therefore we 
determine the soft factor as: 



S(z ,b,( u ,n) 
S(k + ,k ± ,(,LL) 



Re[S 2 (Cu,IJ>)] 



Re 



1 



(2nf 



dz-d 2 be ik+ " 



-ik 1 b 



S(z ,6,C«,/u). 



(36) 



It should be noted that Si and £4 depend on but the soft factor as the ratio of them does not 
depend on £5. With the defined soft factor, the form factors can be factorized as in Eq.(21). They 
take a compact form in the 6-space: 



F V = F A 



iQq 



2v ■ p b^O 



lim / dk + dl + (l) + (k + ,b,C,fJ')S{l + ,b,C u ,n)e(k + +l + )H(k + + l + ,( u ,fj,). (37) 



The limit 6^0 should be taken after the integrations. With the results presented before, we 
determine the perturbative coefficient function H as: 



H(k' 



,Cu,/i) = £jT 



1 + 



2a s (fi) 



3vr 

1, 2k-p 
-I — In — 

2 a(k+r 



2k -p 



2~ Q{k + ) 2 
2k -p 1 




+ Ofa 2 ), 



(38) 



which is free from any soft divergence. All soft singularities are cancelled on a diagram-by-diagram 
basis. The cancellation on a diagram-by-diagram basis is important for extending the factorization 
beyond one-loop level. General arguments for the factorization at any loop can be given by per- 
forming an analysis of relevant reduced diagrams and infrared power-counting. The perturbative 
coefficient function H here does not contain the double log In 2 /j 2 in contrast with that in the 
collinear factorization |14l 1151 116j . instead of ln 2 /x 2 it contains In 2 ( 2 and other log terms. All of 
those log terms need to be resummed if they can be large. 

It should be noted that for the case studied here one may redefine the TMD light-cone wave 
function by including the soft factor as <f>'< = 0+ <g>S, so that the form factors take the form <fi',®H. 
Then our results look similar to those in the collinear factorization. However, it is not clear if the 
same can be done for other processes, because they have not been studied yet. Therefore, we leave 
the soft factor there explicitly. 
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5. Relation between Two Factorizations 



In the section we show that the result of the collinear factorization for the decay can be obtained 
from that of TMD factorization, which is given in the last section. Hence, a simple relation between 
two factorizations is found for the decay. 

In collinear factorization, the transverse momenta of partons are integrated out and the collinear 
light-cone wave function can be defined as|19j: 

M* + ./4 = J ^e ik+Z ~ (01^-^4(00^-^7+75^(00,0)^(0)15^)), (39) 

with the gauge link L n defined with the light-cone vector = (0, 1, 0, 0): 

L n (oo, z) = P exp ^— ig s J dXn ■ G(Xn + z)j . (40) 

By taking the same partonic state as given in Sect. 2., the wave function can be calculated with 
perturbative QCD. The result can be found in ^Uj. With this result and that in the last section, 
one can easily derive the result in the collinear factorization: 

F V = F A = ^L- ( dk + <t> + (k+,v)H c (k + ,v), (41) 
2v ■ p J 

where H c is the perturbative coefficient function and is given by: 




Q(k+) 2 M 2 ( 2 (k + ) 2 C 2 (k + ) 2 

+ 0(a 2 ), (42) 

where the logarithmic terms agree with those in j!41 1151 116j . The difference in constant terms is 
caused by that we used HQET for T^, while full QCD was used to calculate it in [TH ITdI HE] . 
The TMD light-cone wave function has a factorized relation to in b-spacefTU|. It reads: 

/•CO 

4> + (k + ,b,C^)= dq + C B (k + ,q + ,b,(,n)<S> + (q + ,n)+0(b), (43) 
Jo 

where the function Cb can be determined by perturbative theory and is free from any soft diver- 
gence. At leading order of a s the function CB(k + , q + ,b, (, fj,) is 5(k + — q + ). The result of Cb at 
one-loop level can be found in ^U]. It should be noted that from the results in Sect. 2. the TMD 
light-cone wave function <p + (k + , k±, Q at one loop order in the momentum space contains various 
infrared divergences. Some of them are proportional to the tree-level result, i.e., to 5 2 (q±), some 
of them take a form like l/(qj_ + A 2 ). These singularities do not cancel if q± goes to zero. But, 
when we transform 4>+(k + , k±, £) into the 6-space, i.e., when we integrate over k±, some of these 
singularities are cancelled, the remaining singularities are just the same as those in $ + . Therefore 
Cb is free from any soft divergences. The same also happens to the soft factor S with the difference 
that the infrared singularities are completely cancelled, if we transform it into the 6-space, or we 
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integrate over the transverse momentum. The soft factor S in 6-space reads: 



<%+) + ^8tf 



2a, 
' 3vr 



Inf&V) (in | - l) + y + ~ ln 2 (6 2 (P + ) 2 C, 2 ) + 0(b% 



(44) 



with b = be 1 /2. Here S(q + ,b,(u, A*) should be taken as a distribution for q + < P + . The heavy 
quark limit implies P + — » oo. As discussed before and in ^U], we should take finite P + in the 
calculation and take the limit in the final result. The same also applies for Eq.(43), where the 
upper bound of q + should be taken as P + . With these results our factorization formula can be 
re- written as: 

iQ q r p+ 



F V = F A 



2v 



+ 



-Jim lim / dq + <S> + (q + ,ri{HW(q + ,Cu,fi)+H( 1 \q + ,C u ,ri 

■ p b->0 P+^oo Jo L 



P+ 



dk + C ( i ) (k + 1 q + ,C,fi)H^(l + ,C u ,^) 



+ dl+SV(l + ,b,(u,v)H {0) (l + +q + Xu,^ + 0(a 2 s ). 



(45) 



With our results of , <§W and we reproduce H c in Eq.(42). Therefore, the two factorizations 
with fixed orders of perturbative theory are equivalent. 



6. Resummation of Large Logarithms 



In general, one expects that the most important /c + -region of (f) + (k + , kj_, £) for a convolution 
with the wave function like Eq.(37) will be around k + ~ Aqqd- Also the important region of the 
soft factor S(l + ,lj_,(u) is with small l + , i.e., l + ~ Aqcd- This results in that H^\k + + Z + ) will 
contain large single logarithms and large double logarithms and it spoils the perturbative expansion 
of H. Those large logarithms should be resummed for a reliable prediction. 

In the collinear factorization in |14 [ 1151 ITT)] , the resummation can be done by introducing a jet 
factor in the frame work of the soft collinear effective theory [22]. or a jet factor in the full QCD[T%]. 
Similarly, we can also introduce a jet factor in our factorization for the resummation. However, as 
we have seen before, our TMD light-cone wave function and soft factor depend on the parameter 
Q u . This dependence can be used to resum those large logarithms. Before showing this, let us first 
study the evolution of the soft factor. 

The evolution with the renormalization fx and with the parameter C, u reads: 

C 2 



d ~ 
d ~ 

C,u-^rS{k + : b,Cu,n) = 
d 




-^ln<±^ </, + (A ; +,6,C,/x) + 0(a 2 ), (46) 
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where we also include the evolutions of the wave function for completeness. With these equations, 
one can show that the form factors in Eq.(21) or Eq.(37) are independent of as expected. Also, 
their ^-dependence is compensated by the /^-dependence of f(fi) in Eq.(15) so that the matrix 
element (7(6* ,p)\qj fl (l — ^)b\B{v)) does not depend on /x. 

To resum the large log terms, we first take an initial value ( u = (uo m Eq.(37) so that there are 
no large log terms introduced by £m0 in the wave function and the soft factor. Then there will be 
large log terms with ( u q in the coefficient function H, which can be re-expressed with little algebra 
as: 



H(q + , (no,!*) 




In 



2p ■ v 





-) 2 






5vr 2 " 




~ ~6~ 


1 




(47) 



where q + = k + + l + . For small k + and l + there are large log terms in the first line. These terms 
can be resummed by using the C„-evolution Q f H: 



d 

Cu-^—H(q + ,C u ,fj,) 
Solving this equation we have: 

H(q + ,(uO,v) = exp 



2a s 
3vr 



In 



+ ^2 



I 1 



H(q + ,( u ,f,)+0(a 



(48) 




(49) 



All large log terms due to small +-momenta are now resummed in the exponential. To eliminate 
the large log terms in the second line of the above equation, we can set [i = [iu with a large [iu 
so that the ln(2p • v/(ih) is a number of order 1. Then there are large log terms due to large \xh 
in the wave function and the soft factor. With the evolution equations in Eq.(46), we can evolute 
them to lower scales as fi = Finally we have for the form factors: 



F V = F A 



iQ 



^lim / 

• p b^O J 



2v 

. e s• F (fc++^^ 



dk+dl + Mk + , b, k+( u0 , ti )S(l + , b, C u o,^)9(k + + /+) 



k+ + 1+ 



1 + 



2a s {HH) 
3vr 



- In 



2p • v 



5vr2 
6 



(50) 



with 



S F (q + ) 




(51) 



In the above all large logs are resummed in the factor e Sp . The initial value no should be taken 
where perturbative QCD is still applicable. One may take /xq = 1 ~ 2GeV. For ( u0 , with our 
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definitions of the TMD light-cone wave function and the soft factor, we should have £ u q >> 1) 
although the physics here, i.e., the form factors, does not depend on ( u q. However, one should not 
take a too large £ u o to avoid large log terms in the wave function and the soft factor. A detailed 
study of a reasonable choice of CuO and fJ-o is needed when one uses the factorization results for 
phenomenological applications. 

For the resummed results one can also use the relation in Eq.(43) and the result in Eq.(44) 
to express them in term of the usual light-cone wave function $ + , as in the last section. Then 
instead of the integrand <S> in Eq.(45) we have a complicated integrand £g> (g) e Sp . 
Unfortunately, we are unable to calculate the integral analytically. The same also applies to the term 
corresponding to the term in the third line of Eq.(45). Here, we only remind that our resummed 
form factors can be expressed as a convolution of with other functions. 

7. Conclusion and Outlook 

As mentioned in the introduction, there are two approaches for exclusive B-decays. The two 
approaches are not only different in their formulations but also in some predictions in comparison 
with experiment. This leads to controversial discussions, e.g., see I14| 1231 124] . Since two 
approaches are from one fundamental theory-QCD, there must be some relations between them. 
With a consistent definition of TMD light-cone wave functions these relations can be explored and 
predictions for exclusive -B-decays from the two approaches may be unified. For this purpose, a first 
step is to define the TMD light-cone wave function consistently and to obtain relations between the 
TMD light-cone wave function and the usual light-cone wave function in the collinear factorization. 
This has been done in our previous work|lUj. 

In this paper, we have shown that with the consistent definition of the TMD light-cone wave 
function the TMD factorization for the radiative leptonic B-decay can be performed consistently at 
one-loop level. In this factorization, beside the wave function as a nonperturbative object, another 
nonperturbative object, which is the soft factor, must be introduced, so that the perturbative 
coefficient function is free from any soft divergence. The results are given in Eq.(37) and Eq.(38). 
An extension of our factorization beyond one-loop level is possible. 

The TMD light-cone wave function defined in JO] does not only depend on +- and transverse 
components of parton momentum, but also depends on the parameter C, u which regularizes the light- 
cone singularity. This (^-dependence can be used to resum large logarithms in the perturbative 
coefficient function, as we have shown in Sect. 6. 

For the decay studied here, we can show that the result of collinear factorization can be derived 
from that of our TMD factorization. Hence the two factorizations are related to each other. This 
simple relation is to be expected because there is no hadrons, or partons in the final state. In other 
cases, the relation can be expected to be complicated. 

Having shown that TMD factorization works in the simple case, we are ready to explore how 
TMD factorization works in other complicated cases and how it is related to the collinear factor- 
ization in these cases. Works for this are in progress. 
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